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A better understanding of transverse momentum (fc_i_-) dependent quark distributions in a hadron 
is needed to interpret several experimentally observed large angular asymmetries and to clarify the 
fundamental role of gauge links in non-abelian gauge theories. Based on manifestly non-local gauge 
invariant quark operators we introduce process-independent fcx-distributions and study their prop- 
erties in lattice QCD. We find that the longitudinal and transverse momentum dependence approx- 
imately factorizes, in contrast to the behavior of generalized parton distributions. The resulting 
quark fcx-probability densities for the nucleon show characteristic dipole deformations due to cor- 
relations between intrinsic k± and the quark or nucleon spin. Our lattice calculations are based on 
Nf =2+1 mixed action propagators of the LHP collaboration. 



Introduction. — Already 30 years ago, it has been 
noted that intrinsic transverse momentum, fcj_, of par- 
tons gives rise to azimuthal asymmetries in unpolarized 
semi-inclusive deep inelastic scattering (SIDIS), for ex- 
ample e~+p— >e~ +ir+X , nowadays known as the Cahn 
effect [1]. Since then, significant progress has been made 
in understanding intrinsic k± effects and their relation 
to the eikonal phases that quark fields acquire in hadron 
scattering processes due to initial and final state inter- 
actions [2]. The eikonal phases, given by gauge links 
(Wilson lines), turn out to be process-dependent and lead 
to, e.g., the Sivers and Collins asymmetries [3, 4] in po- 
larized S1D1S, which have attracted a lot of attention 
and were already observed in experiments at HERMES, 
COMPASS and Jefferson Lab [5]. Theoretically, these 
can be described in the 
framework of QCD factoriza- 
tion using transverse momen- 
tum dependent parton distri- 
bution functions (tmdPDFs) 
[4, 6], an approach that goes 
beyond the usual collinear 
approximation and operator 
product expansion involving 
(moments of) PDFs. In addi- 
tion to their phcnomcnolog- 
ical importance, tmdPDFs 
provide essential information 
about the internal structure 
of hadrons in the form of 

probability densities in the transverse momentum plane, 
p(x,k±), as illustrated in Fig. 1 [7], where x is the lon- 
gitudinal momentum fraction carried by the quark. 

In this work, we introduce process-independent k±- 
distributions and calculate these in lattice QCD. We il- 
lustrate our results by presenting fc^-densities of quarks 
in the nucleon, with a focus on possible correlations be- 
tween k± and the transverse quark and nucleon spins, 
resulting in deformations from a spherically symmetric 




FIG. 1: Illustration 
of the transverse mo- 
mentum distribution of 
quarks in the proton. 



distribution. It is interesting to compare this approach 
with generalized parton distributions (GPDs) in impact 
parameter (b±-) space [8], which allows one to study the 
spatial distribution of partons in hadrons in form of prob- 
ability densities p(x,b±) [9], Lattice QCD studies of the 
latter revealed characteristic non-spherical shapes of the 
pion and the nucleon in the case of transversely polar- 
ized quarks [10, 11]. We stress, however, that tmdPDFs 
and GPDs provide fundamentally different and comple- 
mentary insight into hadron structure, since they are not 
related by Fourier transformation and k± and b± are not 
conjugate variables. 

To introduce the different tmdPDFs, we first define the 
momentum-space correlators < &r= < i ) r(^: k±; P, S), 
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with nucleon states \P, S) depending on momentum and 
spin, and where the Wilson line U=Uc(ifi), defined by a 
path ordered exponential, ensures gauge invariance of the 
non-local quark operator q(l) . . . q(0). For the vector (un- 
polarized), Ty—j 1 - 1 , axial- vector (polarized), r^=7 M 75, 
and tensor (quark helicity flip), Yl^=ia^ v ^i^, cases, the 
correlators in Eq. 1 can be parametrized by the twist-2 
tmdPDFs [12]: 
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where the distributions f,g,h depend on x and fcj^ and 
A is the nucleon helicity. The light-cone vectors n and 
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n in Eqs. 1,2 are chosen such that n-n—1, P=P + n + 
m 2 N / {2P + )n, and the contraction with n M in Eq. 2 
projects on leading twist for large P + . Formally [7], three 
of the tmdPDFs are directly related to the well-known 
unpolarized, f\{x), polarized, gx(x), and transversity, 
hi{x), PDFs by [fi, gi,hi](x)= f d 2 kj_[fi, gi,hi](x,kj_). 
Further (approximate) relations between tmdPDFs and 
GPDs have been established based on a dynamical mech- 
anism [13], by analogy comparing fej_- and fo^-densities 
[9] , and in the framework of quark models [14] . Such re- 
lations exist in particular for the Sivers, f^~ T , and Boer- 
Mulders, h^, functions, which are naively time-reversal 
odd. Remarkably, the tmdPDFs gyj< and h\ L cannot be 
directly related to any of the GPDs for reasons of time 
reversal symmetry, so that their appearance may be seen 
as genuine sign of intrinsic k± of quarks. In particular, 
they cannot be generated dynamically from coordinate 
(6j_-) space densities by final state interactions. We em- 
phasize that the definition of the correlator in Eq. 1 and 
therefore the tmdPDFs depends on the gauge link path 
C(l, 0). Apart from the precise form of C(l, 0), a further 
issue in the QCD factorization of SIDIS and Drell-Yan 
processes is the appearance of so-called rapidity diver- 
gences, which require regularization or subtraction and 
may be seen in analogy to the Bloch-Nordsieck theorem 
(see [15] and references therein). In the case of SIDIS, 
the path is generically given by C(l, Q)=[l,l + oon, oon, 0], 
which is not invariant under time reversal, and hence 
the appearance of the Sivers and Boer-Mulders functions. 
There exist two natural choices for process-independent 
fc_i_-distributions, employing either a direct Wilson line 
C(l,0)=[l,0], or an average over all possible gauge links 
between and /. We chose the first possibility for this 
pioneering study. Our results should not be directly com- 
pared to phenomenological studies of SIDIS [16] or the 
Drell-Yan process. In particular, the Sivers and Boer- 
Mulders functions vanish for direct Wilson lines. To de- 
scribe the corresponding asymmetries, we will have to 
use specific gauge link structures, which is the subject 
of ongoing studies [17, 18]. Lattice QCD gives access to 
hadron matrix elements <I>r(Z;P, S) (cf. Eq. 1), which 
can be parametrized by complex- valued invariant ampli- 
tudes Ai(l 2 ,l-P). The choice of a straight Wilson line, 
together with hermiticity, parity and time-reversal sym- 
metry leaves us with 8 independent amplitudes: 
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where \jni\=jiv — v[i. Detailed comparison with Eqs. 2 
reveals that the tmdPDFs are given by certain Fourier- 
integrals of the amplitudes, where the Fourier conjugate 
variables are l-P and x, and lj_ and k±. Such relations 



exist between fi<->A 2 , gx<->A 6>7 , hx<->A 9m , g XT ^A 7 , 
h^ L ^Ai n, and h^~ T <->An, while ^3,8 give only con- 
tributions to higher-twist tmdPDFs. 

Lattice QCD results. — In our lattice study, we use a 
discretized, Euclidean version of the non-local operator 
Op(l)=q(l) TUuo] <?(0). Limiting ourselves to purely spa- 
tial quark separations, Z°=0, the operator resides on a 
plane of equal Euclidean time. A product of connected 
link variables is used to represent the straight Wilson 
line U[ifl]- For oblique angles, we approximate a straight 
line by a step-like (zig-zag) path. According to the con- 
tinuum analysis of Ref. [19], the renormalized operator 
can be written as Op n (/)=Z _1 exp(— Sm yj—l 2 ) Or(l), 
involving two renormalization constants, Z^ 1 and Sm, 
the latter being associated with a potential power diver- 
gence of the Wilson line. The renormalization of Wil- 
son lines in a cut-off regularized theory, e.g. on the lat- 
tice, and the appearance of power divergences of the form 
I /a has been studied extensively in the context of heavy 
quark effective theory, see, e.g., [20]. We utilize the fact 
that Sm also appears in the renormalization of the static 
quark potential: V rcn (r)—V(r) + 25m. At large r, the 
static quark potential is well approximated by the string 
potential V s \™ ng (r)=ar - n/l2r + C* rcn . To fix Sm we 
follow Refs. [21] and match V ren (r) to the string po- 
tential at r = 1.5r «0.7fm, setting C rcn =0. Note that 
fc_i_-distributions are sensitive to the renormalization con- 
dition implied by the choice of C rcn . In our numeri- 
cal computations, we employ MILC gauge configurations 
[22] based on the AsqTad improved staggered quark ac- 
tion with two light and one heavier (strange) quark fla- 
vors. For our exploratory calculations, we have chosen 
an ensemble at a pion mass of «500MeV and a lattice 
spacing of a=0.124fm. From an analysis of V(r) from 
Wilson loops on the HYP smeared gauge configurations, 
we obtain aSm= — 0.155(5) sta t [17, 18]. To determine 
the amplitudes Ai, we take advantage of the methods 
and techniques used in Ref. [23] for the calculation of 
GPDs. As in Ref. [23], we only include contributions 
from quark line connected diagrams. A major innova- 
tion in our analysis is that the lattice correlation (three- 
point) functions are evaluated for manifestly non-local 
operators Or (I), hence requiring a parametrization ac- 
cording to Eq. 3. They are numerically evaluated using 
the CHROMA library [24] for domain wall quark propa- 
gators and sequential propagators previously calculated 
by the LHP collaboration, with the valence quark mass 
tuned to match the staggered sea quarks [23]. The se- 
quential propagators are available for two lattice nuclcon 
momenta P=(0,0,0) and P=(— 1,0,0), the latter corre- 
sponding to «500MeV in physical units. We consistently 
use the nucleon mass of mjv=l. 291(23) GeV obtained for 
this ensemble in the extraction of the amplitudes and 
tmdPDFs from Eqns. (2,3). All errors quoted are statis- 
tical (see Refs. [17, 18] for a more comprehensive treat- 
ment of uncertainties). Setting l°—0, we are restricted to 
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FIG. 2: Real parts of the amplitudes Ai(l 2 , l-P=0), for 
m^~500MeV. The solid lines and error bands are Gaussian 
fits to the renormalized lattice data above %/ — Z 2 =0.25 fm. 





c 


2/a (GeV) 




9 fll 59(86) 


_ f(0,0) 
— Jl,u 


0.3741(72) 


At 


1 0199(90) 


_ f(0,0) 
— Jl,d 


0.3839(78) 


At 


-fl 990(^5) 


- -o (0 ' 0) 


0.311(11) 


At 


n 9Qi (i 9) 


- -o (0 ' 0) 


0.363(18) 


A u 


CHI 


_ 1,(0,0) 
— a l,u 


0.3184(90) 


M 


-0 254(16') 


"l,d 


0.327(15) 


Ay 


-0.1055(66) 


- -o ( ° !l) 


0.328(14) 


Af 


0.0235(38) 


- -o (0>1) 


0.346(36) 


ifij 


-0.0931(73) 


— a lL,u 


0.340(14) 


Mo 


0.0130(40) 


n lL,d 


0.301(48) 



TABLE I: Parameters obtained from a fit of the amplitudes 
Ai{—l\,Q) to Gaussian functions, for m^~500MeV. 



scription at short distances I 2 , which translates into a 
smooth cutoff of the resulting fc^-distributions at large 
k\. Within this context, we define fe^-moments 



the region given by Z 2 < and \l-P\ < V~l 2 \P\- Since 
only finite hadron momenta P can be accessed on the lat- 
tice, this precludes us from evaluating the full x- and /in- 
dependence directly. An important result, demonstrated 
in detail by our lattice data in [17, 18] in the kinematic 
region accessible to us, is factorization of the x and k± 
(l-P and I 2 ) dependence to a good approximation. Fac- 
torization is a standard assumption for which we have 
now provided quantitative support, and for this reason 
we will focus in the following on the I dependence and 
present results for the lowest x-moments corresponding 
to l-P = 0, e.g., 
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We note that the ^-integrated distributions of type 
and h^ ' correspond to the difference, while distributions 
of type f/ ' correspond to the sum of quark- and an- 
tiquark distributions, respectively. Hermitian conjuga- 
tion shows that lmAi(l 2 , l-P— 0)—0. Representative re- 
sults for ReA 2! 7 i io(— ^1, 0) are displayed in Fig. 2. Sim- 
ilar results not showii^in this work have been obtained 
for the amplitudes AQ i g m> ±i. We have parametrized 
our renormalized lattice data using a Gaussian ansatz 
2^«, (— / 2 L) 0)=c"' exp(— ix/cj*' ). Corresponding fits 
to the lattice data are represented by the error bands in 
Fig. 2. To avoid contamination by lattice cut-off effects 
at small l^a,we exclude data points for V—l 2 < 0.25 fm 
from the fit. The Gaussians smoothly bridge the ex- 
cluded region at small I 2 . Thus we do not resolve po- 
tential divergences of the amplitudes Ai(l 2 ,0) in the 
continuum limit at I 2 — >0. In that sense, the Gaussian 
parametrization effectively acts as a regularization pre- 
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Since f[ (fcj_) ^ s the unpolarized density of quarks mi- 
nus antiquarks, we can fix Z by demanding that the 
total number of u- minus c?-quarks in the proton is 

fi°u-d= c 2~ d = !> g ivin S 2=1.056(14). The Gaus- 
sian fit parameters are listed in Table I. Note that the 
width a is very sensitive to our renormalization condition 
C ron = 0. For A 2fi , 0m , the inverse width 2/a^(k 2 L ) 1 / 2 
has an interpretation as the root mean square trans- 
verse momentum of the respective distribution /} ' , , 



hf\ The fact that gr^_ d =1.209(36) comes out close 
to the physical gA—^- 2695(29) and agrees within errors 
with the direct calculation in [25] represents a remark- 
able validation of our approach. Interestingly, we find 
that gj^ 1 '^ — hfj^ ' 1 ^, which supports corresponding re- 
sults obtained in quark model calculations [26]. The 
fej_-densities of quarks in the nucleon for longitudinally 
(L) and transversely polarized quarks (T) can be ob- 
tained from pi J —p(x,k±;\,S±)—n-(^v + A ( f> J 4)/2 and 
pr=p(x, k±; A, s±)=(n-<&v — n^Sj$tf')/2, respectively, 
and appear in form of a multipole-expansion [9] 
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with monopole terms oc fi,gi,hi, dipole structures ex 
/lTjffiT) hi , h± L , and a quadrupole term oc hj~ T . The 
terms in square brackets proportional to the T-odd Sivers 
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FIG. 3: Quark densities in the fc^-plane, for m^~500MeV. 
(a) pL for u-quarks and A = 1, S± = (1,0), (b) the same for 
d-quarks, (c) pr for it-quarks and A = 1, Sj_ = (1,0), (d) the 
same for d-quarks. The error bands show the density profile 

(scale not shown). 



and Boer-Mulders functions are absent for straight Wil- 
son lines. We stress that the dipole-correlations of the 
types oc Xk±-S±giT and oc Ak±-s±h^ L in Eqs. 4 are 
a characteristic feature of intrinsic transverse momen- 
tum, and that analogous terms in the case of spatial 
(impact parameter) densities are forbidden by time re- 
versal symmetry [9]. Based on our results for the k±- 
distributions, we show in the upper part of Fig. 3 the 
lowest i-momcnt of the fc^-density pl of longitudinally 
polarized u- and d-quarks in a transversely polarized nu- 
cleon, with A= + 1 and iSj_=(l,0). The densities feature 
significant dipole deformations along the direction of the 
nucleon spin Sj_ , due to non- vanishing average transverse 
momentum shifts (k x )=mpjgf > j, 1 ' > / f[°'°\ We find sizeable 
shifts of (fc a .)=67(5) MeV for u-, and (fe x )= - 30(5) MeV 
for d-quarks. An analogous observation is made for trans- 
versely polarized quarks in a longitudinally polarized nu- 
cleon, as displayed in the lower part of Fig. 3 for A= + 1 
and sj_=(l) 0)- Here we find shifts of similar magnitude 
but opposite sign: {k x )=m^h^' 1 ^ / '°^= — 60(5) MeV 
for u-, and (k x )=16(5) MeV for d-quarks. We note that 
these intrinsic fcj_-shifts are orthogonal to the dipole de- 
formations oc SjEjibi of densities in impact parameter 
space observed in [10]. 

Conclusions. — We studied tmdPDFs of the nucleon 
for the first time using lattice QCD. We found that 
the x and the k± dependence approximately factorizes 
and that quarks in the nucleon carry significant intrinsic 



transverse momentum. The fc_L-densities for polarized 
quarks in a polarized nucleon show characteristic dipole 
deformations along the transverse spin vectors. This 
work represents a first milestone in the non-perturbative 
calculation of fcj_-distributions. Ongoing studies includ- 
ing non-straight Wilson lines promise in particular to give 
access to the highly interesting T-odd Sivers effect. 
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